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Abstract 

It is shown that, in a reasonable approximation, the quantum state of p-bosons in a bi-partite 
square two-dimensional optical lattice is governed by the nonlinear boson model describing tun- 
neling of boson pairs between two orthogonal degenerate quasi momenta on the edge of the first 
Brillouin zone. The interplay between the lattice anisotropy and the atomic interactions leads 
to the second-order phase transition between the number-squeezed and coherent phase states of 
the p-bosons. In the isotropic case of the recent experiment. Nature Physicis 7, 147 (2011), the 
p-bosons are in the coherent phase state, where the relative global phase between the two quasi 
momenta is defined only up to mod(7r): (j) = ib7r/2. The quantum phase diagram of the nonlinear 
boson model is given. 
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Cold atoms and Bose-Einstein condensates in optical lattices provide a versatile tool for 
exploration of the quantum phenomena of condensed matter physics on one hand, and, 
on the other hand, a way for creation of novel types of order in cold atomic gases [1]. 
Two remarkable recent achievements in this direction are the experimentally demonstrated 
novel types of atomic supcrfluids in the P- [2] and F-bands [3] of the bi-partitc square 
two-dimensional optical lattice. The bi-partite optical lattice having a checkerboard set 
of deep and shallow wells (i.e. made of double- wells) , used in Refs. [2, 3] has a large 
coherence time in the higher bands, several orders of magnitude larger than the typical 
nearest-neighbor tunneling time [4]. The order parameter of these superfluids is complex, 
in contrast to the conventional Bose-Einstein condensates having real order parameter in 
accord with Feynman's no-nodc theorem for the ground state of a system of interacting 
bosons [5, 6]. The p-bosons, for instance, are confined to the second Bloch band for a 
sufficiently shallow lattice amplitude, Vq < 2.2Eji, where is the recoil energy [4, 7]. In 
Ref. [2] Vq ~ 1.55Eji, however, a particular experimental technique was used which results 
in population of other Bloch bands. Nevertheless, the main results on the cross-dimensional 
coherence are obtained for the parameter values where the second band is by far the largest 
populated. 

The purpose of this work is to show that, in the reasonable approximation, the quantum 
state of the p-bosons in the square bi-partite optical lattice is governed by the modified 
nonlinear boson model, which was already used before in the context of cold atoms tun- 
neling between the high-symmetry points of the Brillouin zone [8-10]. However, there is 
an important difference: in the p-boson case there is a lattice asymmetry parameter which 
provides for the phase transition at the bottom of the energy spectrum, additionally to that 
at the top of the spectrum, studied before in Ref. [9]. The focus is on the quantum features 
of the p-boson superfluid, as different from Ref. [11] where a mean-field Gross-Pitaevskii 
approach was employed and the region of the complex order parameter was found. 

The nonlinear boson model derived below follows just from two basic conditions: the 
existence of two quasi degenerate energy states coupled by the boson pair exchange (tun- 
neling) when the single-particle exchange is forbidden. Thus it applies to other contexts as 
well (see also Ref. [10]). For instance, it is equivalent to the nonlinear part of the so-called 
fundamental Hamiltonian (in the Wannier basis) , describing the local two- flavor collisions in 
the first excited band of a two-dimensional single- well optical lattice [12] . Moreover, in the 
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case of the optical lattice consisting of the one-dimensional double- wells [13] the many-body 
Hamiltonian can be cast as a set of linearly-coupled nonlinear boson models. Taking this 
into account, we consider the quantum features of the derived nonlinear boson model in the 
most general setting and using its natural parameters, besides analyzing the experimental 
setting of Ref . [2] . 

Consider the bi-partite square two-dimensional optical lattice of Ref. [2], which can be 
cast as follows (after dropping an inessential constant term) 

V = — Voexp^ — (2^2,0 cos(2A;a;) -|- 2^0,2 cos(2A;y) 

+2Re{T;i,ie''=(^+^) + T;_i,ie^*^(^-^)}) , (1) 

where the experimental values of the parameters read Vq — Vo/4 — 1.55Eji in terms of the 
recoil energy — with k — 2n / X, X — 1064 nm and — 71 //m being the oscillator 
length of the transverse trap. The dimensionless Fourier amplitudes of the lattice are 

"^2,0 — cos a, vi^i — rje [e*^ -|- e~*^ cos a] , 

Vo^2 = e, ^^-1,1 = V [e'^ cos a + e^e~'^] , (2) 

see Fig. 1. The experimental parameters are rj 0.95 and e ^ 0.81. 

For aiso — arccose ~ 7r/5 and arbitrary values of the other parameters we have vi^i — 
hence, the lattice satisfies the symmetry V{—x,y,z) — V{x,y,z). For a — a^gg the 
band energies with the Bloch indices Ki = (|, |) and K2 = (— |, |) (sec Fig. 1(a)) become 
equal, since the Bloch functions satisfy LpiQ-^{—x,y,z) = LpiQ,^{x,y, z), due to the boundary 
conditions and symmetry of V{x, y, z). Hence, the points are the high-symmetry points 
of the symmetric lattice with (see also Ref. [8]). 

As is found in Ref. [11] the observed cross-dimensional coherence [2] is the joint effect 
of the atomic interactions and the lattice potential. Indeed, let us estimate the interaction 
energy and its characteristic time scale in the p-boson experiment. The interaction energy 
can be estimated as Eint ~ where g = ^^^"^ is the interaction cocfScicnt proportional 
to the s-wave scattering length a^, iV is the number of atoms and Q is the effective volume 
of the condensate. Setting fl — V^a^L^A^ (the coefficient is due to ground state of the 
transverse trap, see below), where L is the number of the lattice sites along each of the 
two directions di_2 in the plane x = (x, y), we get E^nt ~ ^^^R- ^^Rb and other 
experimental values of Ref. [2] , with N ~ 10^ and L ~ 10 (the estimated sample size of Ref. 
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FIG. 1: (Color online) Panel (a): The Fourier spectrum of the 2D lattice (1), where the Fourier 
amplitudes are shown by the (blue) dot markers and pointed by the reciprocal lattice vectors, the 
vectors Qi = {k, k) and Q2 = {—k, k) give the reciprocal lattice periods, the points Ki^2, shown by 
two (red) dot markers, are two nonequivalent minima of the second Bloch band. The square area 
is the first Brillouin zone. Panel (b): The bi-partite square lattice in the real space, the vectors 
di = (^, and d2 = (— f , f ) form the basis of the lattice periods. 

[2] divided by the lattice cell size), we get Eint/Eji ~ 0.1. Moreover, the interaction time 
scale, defined as tint = fi/Eint is on the order of the typical experimental times, indeed, we 
have tint 

~ ^'m^sn' ~ 10ms (compare with Fig. 2 of Ref. [2]). 
Taking into account the above estimate, one can assume that the atoms are confined to 
the second Bloch band of the lattice and expand the boson field operator over the band- 
limited Bloch basis. The Bloch waves are defined as V5k(x) = -^e*''''uk(x), Mk(x+dj) = Uk(x), 
j = 1,2, where the periodic Bloch functions Mk(x) are chosen to be normalized on the 2D 
lattice cell z/q = |di x of area A^/2, i.e. 

I d2x|«k(x)P = l. 

The band-limited expansion reads 

^(x,2;) = te(x)<l'o(^), (3) 
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where the summation is over the Bloch indices inside the first Brillouin zone k e 
{fQi + fQa; Kj = -L/2,...,L/2} (see, Fig. 1(a)). Here %{z) is ground state 
$0(2) ~ 7r~^/^aJ^^^e~^^/(^"^) of the transverse trap. Inserting this expression into the stan- 
dard Bose-Hubbard Hamiltonian for the lattice potential (1) and using the Poisson summa- 
tion formula, 

L 

^ gik(^idi+^.d.) ^ ^2^^^^ mod(Q), 

we obtain 
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E ^(k)^'^k&k + 1^ E X(ki,k2|k3,k4)6tk,6t^,6k3&k4, 

keBZ Ak=0 



(4) 



where -E'(k) is the Bloch energy of the second band, Ak = ki -|- k2 — — k4, the condition 
Ak = is understood mod(Q) and 



X(ki,k2|k3,k4) = i/Q j d^xuki<2iik3^ik4 



(5) 



is the dimensionless coefficient which depends solely on the lattice geometry 

Since the points Ki 2, the energy minima of the second band, are lying on the edge of 
the Brillouin zone (Fig. 1(a)), the Bloch functions </?Ki,2(x) are real. Moreover Vk-E(k) = 
and, hence, Vk</'Ki,2(x) = 0. As the result, the expansion over k in Eq. (4) in some small 
neighborhoods about these points starts only with the second-order term oc (k — Ki^2)^ 
[14]. On the other hand, one can verify that the experimental width of the Bragg peaks 
about the band minima Ki^2 is too narrow to give a significant second-order correction, i.e. 
(k — Ki^2)^/^^ ~ 0.06 (see. Fig. 3 of Ref. [2]). Therefore, we can discard the spectral width 
of the Bragg peaks and keep in Eq. (3) only the two-mode expansion of the boson field 
operator (a similar expansion over the two nonlinear modes was also used in Ref. [11]) 



*(X,^) W [6l(^Ki(x) + 62¥'K2(x)] $o(^)- 



(6) 



It is important to note that, since the summation in the nonlinear term of Eq. (4) is 
conditioned by Ak=0 mod(Q), all terms with with either three Ki and one K2, or vice 
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versa are zero (i.e. bosons tunnel between the minima by pairs [8]). Thus, only the following 
geometric parameters are nonzero: 



Xn ^ ^0 / d^^l^K^ (^) ^ i = 1> 2, 

Xl2 = 2^0 j d^x|liKi(x)liK2(x)|^. (7) 



As a consequence, one obtains from Eq. (4) the two-mode Hamiltonian of the nonlinear 
boson model [8-10] except for the term proportional to the population imbalance due to the 
lattice asymmetry: 

£1—82 U 
H = — - — (ni - n2) + ^{ni(ni - 1) + 712(712 - 1) 

+A(47ii7i2 + + {bW)}- (8) 

We have denoted rij = 6^6;. The parameters of Hamiltonian (8) are as follows. The energies 
of the two symmetric points read 

where i^i 2 = -E(Ki 2) is the respective Bloch energy, = 711+712, U is the average interaction 
parameter per particle, 

t/=|^(Xll + X22), (10) 

and A is a pure geometric parameter defined as 

A = 2Xl2/(Xll+X22). (11) 

Note that at the symmetric point a — aiso we have (7 = 0, hence £1 — £2- We have just two 
independent parameters (7, A), where 7 is defined as 

£l-£2 E1-E2 _Xll-X22 

Here we note that any 2D lattice which for some set of parameters possesses two non- 
equivalent points lying on the edge of the Brillouin zone and having equal Bloch energies 
can lead, under similar conditions, to the same model Hamiltonian (8). 

The parameters A and cr, < A < 1 and — 1 < cr < 1, are independent of the interaction 
strength g and are functions only of the lattice shape. For the experimental lattice (1) their 
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dependence on a and 6 can be determined by numerically solving the 2D eigenvalue problem 
for Bloch energies, the result is given in Fig. 2. Except for the semicircle shaped plateau, 
both parameters vary significantly with variation of the lattice potential. Specifically, for 
the experimental value 6 = O.SStt the parameters A, a and the Bloch energy difference are 
given in Fig. 3. 








FIG. 2: (Color online) The numerically computed lattice parameters A (panel (a)) and a (panel(b)) 
for the experimental 2D lattice (1) of Ref. [2]. Here Vq = 1.55Er, r] = 0.95, e = 0.81. (The high 
accuracy Fourier pseudospectral method [15] was used.) 

The interaction energy parameter UN was already estimated above, i.e. UN/Ept ~ 
ttsN / (ozL^) ~ 0.1 for the experimental values of Ref. [2], the bandgap is on the order of the 
lattice amplitude Vq ~ Er, o ~ 0.2, see Fig. 2(b), whereas the energy degeneracy is at most 
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FIG. 3: (Color online) The energy difference, (a), and the lattice parameters A (the left y axis, 
solid line) and a (the right y axis, dashed line), (b). Here 6 = O.SStt and the other parameters are 
as in Fig. 2. 



~ OAEji, see Fig. 3(a). We conclude that in the experiment of Ref. [2] 7 can reach order 
~ 1. 

The ground state of the model Hamiltonian (8) for 7 7^ can be one of the two types 
of states: either the coherent phase state (with a definite relative phase between the two 
modes 61,2, = ±7r/2) or the atom number squeezed (Bogoliubov) state. These two types 
of the ground state are connected by the second-order quantum phase transition on the 
borderlines A = l-f I7I in the plane (7, A) (see also Fig. 4 below). There are, in fact, exactly 
two phase transitions. One is at the bottom of the quantum energy spectrum and occurs 
for the asymmetry parameter 7 7^ (it corresponds to the relative phase ±7r/2, see below). 
The other one is at the top of the spectrum (and corresponds to the zero relative phase). 
For 7 = the phase transition at the top of the spectrum was studied before [9, 10]. 

Consider first the number-squeezed states, which appear for the large population imbal- 
ance between the points Ki 2 and have a squeezed variance of the population imbalance 
(see also Refs. [9, 10]). For instance, suppose that rii ^ ^2 (i.e. rii ^ A^) and denote 
the respective class of states by Bi. Following Bogoluibov's approach, one can replace 
bi — )■ \/N — n2e**, where $ is an inessential random phase, and expand the Hamiltonian 
(8) in orders of 62 and h\. Keeping the second-order terms only we get the local quadratic 
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Hamiltonian in the form H ^^^Hi with 

^1 = (1 + 7) + 2(2A - 1 + 7)^ + ^[{e^^blr + {e-'^b,n (13) 
The Hamiltonian (13) is diagonahzable by the Bogohubov transformation 

62 = e^*[cosh/3a - sinh/3a^], tanh(2/3) = — ■ (14) 

2A — 1 — 7 

where (3 is the squeezing parameter. We have 

2 / . 1 



^1 = 1 + 7 + ^ 



a^a + - coth(2/3) 



(15) 



_sinh(2/3) 

For ^2 ^ rii the number-squeezed states S2 are described by similar quadratic Hamiltonian 
H2 obtained by replacing 62 with bi in Eqs. (13)-(14) as well as inverting the sign at 7. 

The existence diagram of the number-squeezed states B12 is shown in Fig. 4(a), their 
existence is equivalent to existence of the Bogoliubov transformation (14). The states Bi^2 
are thermodynamically stable for positive effective mass in Eq. (15), i.e. when 2A — 1=1=7 > 0) 
which condition is satisfied only in the regions A > 1 -|- 7 and A > 1 — 7, respectively for Bi 
and B2. The thermodynamically stable Bi^2 states are shown in Fig. 4(d). Note that the 
number-squeezed states have undefined relative phase (f) — arg(((6|)^6|))/2 (this is reflected 
also in arbitrariness of see also the discussion of the quantum phase below). 

Hamiltonian (8) also admits the phase states possessing definite values (i.e. with small 
variance) of the phase and the population imbalance. These states will be called coherent. 
The existence diagram of the coherent states can be found by approximating the Hamiltonian 
by a quantum oscillator problem in the Fock space [9] . For N ^ 1 the coherent states are 
essentially semiclassical in the sense of Ref. [16]. Thus, most of their properties can be 
studied by replacing the boson operators by scalar amphtudes: bi — >■ -y/iV72(T+~C)e~*'^/^ 



and 62 ~^ \/N/2{\ — Qe^'^l'^ and considering the resulting classical model (save for the 
factor ^) 

^ci = ^ + 7C + ^ [C' + A(l - C')(2 + cos(20))] . (16) 

The stable stationary points of the classical Hamiltonian "Hci correspond to the phase states 
of the quantum model. There are two stationary points: (20^ = 0, = 3x1:1) ^"^^ (206 = 
TT, Cfe = ~i^) s-iid they correspond, respectively, the coherent phase states at the top (Cq) 
and at the bottom (C^r) of the quantum energy spectrum (this is clear from their energies). 



9 



The direct approach to study the coherent states is based on the discrete WKB in the 
Fock space, with the effective Planck constant h — 2/N [9]. One first factors out the classical 
phase 06 1 and then expands the Hamiltonian (8) about the classical stationary point (b^t (see 
also Ref. [17]). Representing the Fock-spacc "wave function" ip{C) = {n,N — n\ilj) (here 
n = rii) with ^ = 2n/N — 1 as ip = e^'^^^^ipoiC,) and defining the canonical with C, momentum 
as p = —ihd(^ we get 

(n, N - = e^'^^l^'^H^UCl (17) 

with a local Hamiltonian of a quantum oscillator (the discarded terms start with ~ 
(C ~ Qb,tf')- The Hamiltonian about Cfe (for the phase 20^ = tt) reads 

while that about the point Ct (for 20^ = 0) can be obtained by replacing A by 3A in the 
first two terms in Eq. (18) and inverting the sign at due to the negative effective mass 
M^l = —K/2{l — Q^) cos{2(f)b,t)- The existence and stability analysis is straightforward from 
this point. First of all, the coherent states Co, i.e. with the classical phase satisfying 20^ = 0, 
are thermodynamically unstable due to the negative effective mass, while the states are 
thermodynamically stable where they exist. The existence diagram of the coherent states is 
given in Figs. 4(b) and (c). Numerical simulations confirm that the Gaussian width of the 
oscillator "wave- function" ■0(C) reasonably approximates the width of the coherent states in 
the Fock space. 

By considering the characteristic energies (up to ~ 1/A^) in terms of UN'^/2 of all the 
above classes of states, i.e. E{Bi 2) = 1 ± 7, E{Ct,) = ^2(1-/^ ^'^'^ E{Cq) = ^2(1-3 A) ' 
obtains the state diagram of the model (8), see Fig. 4(d). Depending on the values of 7 
and A the ground state is either the coherent state C^r or one of the squeezed states, Bi or 
B2. The phase transition borderhne is A = 1 + |7|. Figs. 4(a) and (b) demonstrate that a 
similar phase transition occurs at the top of the energy spectrum on the border line given 
by A = (1 + |7|)/3. It was the subject of Rcfs. [9, 10]. 

Let us now consider the state diagram versus the experimental parameter a. To compare 
the result also to the mean-field diagram of Ref. [11] (see Fig. 5) one has to identify 
the same interaction parameter (the product of the g and the density in Ref. [11]). The 
quantity gN/fl can serve as an analog, though one has to remember that we have discarded 
the atoms of the condensate not represented by the Bragg peaks at the two points Ki 2, 
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FIG. 4: (Color online) (a) The existence diagram of the number-squeezed (Bogoliubov) states Bi^2 
of Hamiltonian (8). (b) The coherent states of zero relative phase, corresponding to the top of 
the quantum energy spectrum (Co), (c) The coherent states with the relative phase (p = ±7r/2, 
corresponding to the bottom of the quantum energy spectrum (Ct^). (d) The ground state diagram 
of the model Hamiltonian (8). 



thus the resulting approximate value of gN/Q will be smaller than the actual value and 
the comparison can be only qualitative. The expressions for the borderlines A = 1 ± 7 of 
the state diagram Fig. 4(d) can be rewritten using Eqs. (9), (10) and (12) as to give the 
interaction parameter gN/Q. We obtain: 

N E2-E1 

9n=— ^^0' 19 

Q Xu - X12 

N E1-E2 
Q X22 - X12 

The results are presented in Fig. 5, where the energy is given in the recoil energy units. 
Qualitatively we have similar diagram to that of Ref. [11], though the corresponding quan- 
titative value of the interaction parameter gN/Q is significantly smaller (though the density 
parameters are not identical, as mentioned above, the difference is still significant). We 
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note, however, that the values of the interaction parameter in Fig. 5 do correspond to the 
estimated value gN/Q ~ f/A^ ~ 0.1Ef> which accounts, for instance, for the Bragg peak for- 
mation times. This estimate was used to validate the expansion (6) over the Bloch modes, 
which was then used in the nonlinear part of the many-body boson Hamiltonian to produce 
the model Hamiltonian (8). For this very reason only the lower part of the figure around 
the critical aiso belongs to the validity region of the approximation. Finally, an analog of 
the relative populations of the two modes is the semiclassical imbalance (b (defined only for 
the coherent states). It can be cast as 

X22 - Xii + 2(^2 - E,)n/{gN) 



7 

1 - A 



Xii + X22 - 2x 



(21) 
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see Fig. 5(b). 
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FIG. 5: (Color online) (a) The phase diagram of the model Hamiltonian (8) in terms of the effective 
interaction parameter in units of the recoil energy vs. the angle a. Here 9 = O.SSvr and the 
rest of the parameters as in Fig. 2. (b) The typical average relative population imbalance between 
the two points Ki^25 which correspond to the semiclassical between two values of a (the dotted 
lines in (a) and (b)). 



Finally, let us make some comments on the relative phase (f). Why the phase 20 appears 
in the classical Hamiltonian Hd (16) is clear: the bosons tunnel by pairs, which is reflected 
in the splitting of the even and odd subspaces of the Fock space, with the respective basis 
states |2s, N — 2s) and |2s — 1, — 2s + 1) [9, 10]. Since the state of the system is always 
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expanded over the states differing by an even number of bosons, it is impossible to define tlie 
pliase 0, but only the 20: 20 = arg(((6|)^6|)). Hence 20 and not appears in the exponent 
factor in Eq. (17): exp{i(p(/h} = exp{2i0(ni — 712)}. The splitting of the Fock space into 
two subspaccs also leads to the double degeneracy of the coherent states (quasi-degeneracy 
to be precise: the terms of order 1/A^ are neglected), since the same approximate "wave- 
function" in the Fock space ip{C) describes not one but two states, one of each subspace: 
C2s = (2s, - 2s\ijj) and C2S-1 = (2s - 1,A^ - 2s + 1|V') with the discrete sets Ci e 
{(2s -1)/N- 1} and C2 e {2s/A^ - 1}. 

The mean-field approach, in contrast, produces a definite relative phase, see Ref. [11], 
where two equivalent order parameters of the nonlinear Gross-Pitaevskii equation arc pos- 
sible for the description of the same experiment with the phase cither ±7r/2, due to the 
broken superposition principle by the nonlinearity. However, the full many-body quantum 
Hamiltonian permits superposition of the eigenstates of the same energy. The resolution 
of this seemingly paradoxical situation is similar to the case of the random phase in the 
double-sht experiment with the Bose- Einstein condensate, see Ref. [18]. Indeed, since the 
atoms are detected one by one coherently from both modes 61,2, when the lattice is released, 
the atom detections probe the quantity {b\b2) spontaneously projecting, as the detection 
process proceeds, on one of the two possible phases 0^ = ±7r/2 of C,r- 

In conclusion, we have shown that the experiment of Ref. [2] is describable by the 
quantum model (8) and that there is the quantum phase transition of the second order 
between the atom number-squeezed states and the coherent phase states of the p-bosons. 
The results indicate that in the recent experiment [2] a phase transition of the second 
order was observed, where the isotropic experimental state observed for the symmetric point 
a — arccose (and hence, for 7 = 0) must be the coherent Ctt state of the relative phase 

20 = TT. 
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